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Abstract 

Lipid-bilayers are the fundamental constituents of the walls of most living 
cells and lipid vesicles, giving them shape and compartment. The formation 
and growing of pores in a lipid bilayer have attracted considerable attention 
from an energetic point of view in recent years. Such pores permit targeted 
delivery of drugs and genes to the cell, and regulate the concentration of 
various molecules within the cell. The formation of such pores is caused 
by various reasons such as changes in cell environment, mechanical stress or 
thermal fluctuations. Understanding the energy and elastic behaviour of a 
lipid-bilayer edge is crucial for controlling the formation and growth of such 
pores. In the present work, the interactions in the molecular level are used 
to obtain the free energy of the edge of an open lipid bilayer. The resulted 
free-energy density includes terms associated with flexural and torsional en¬ 
ergies of the edge, in addition to a line-tension contribution. The line tension, 
clastic moduli, and spontaneous normal and geodesic curvatures of the edge 
are obtained as functions of molecular distribution, molecular dimensions, 
cutoff distance, and the interaction strength. These parameters are further 
analyzed by implementing a soft-core interaction potential in the microphys¬ 
ical model. The dependence of the elastic free-energy of the edge to the size 
of the pore is reinvestigated through an illustrative example, and the results 
are found to be in agreement with the previous observations. 

Keywords: Free energy, elasticity, open lipid bilayer, differential geometry, 
boundary curve of a surface, molecular interactions. 


* Corresponding author 

Email address: meisam. asgariQmail .mcgill. ca (Meisam Asgari) 


Preprint submitted to International Journal of Non-linear Mechanics 


April 7, 2015 



1. Introduction 


A phospholipid molecule consists of a hydrophilic head and two hydrophobic 
fatty-acid tails [Tj. When suspended in an aqueous solution at sufficient 
concentrations, phospholipid molecules self-assemble into structures such as 
lipid bilayers, in order to shield the tail groups from the solvent M. Lipid 
bilayers are the main constituents of cell membrane in most living organisms, 
as well as model membranes such as liposomes [4], They provide the cell 
and its substructures with compartment and shape, and further, function as 
barriers for water-soluble molecules such as water, ions, and proteins PEE]- 
Lipid bilayers are composed of two adjacent leaflets of phospholipid molecules 
oriented transversely and set tail-to-tail. 

Forming of open edges in lipid membranes results in exposer of the tail 
groups at the edge to water [4], which is energetically unfavourable. As a 
result, phospholipid molecules rapidly rearrange around the exposed edge, 
forming a semicylindrical rim along it. This rearrangement is the source of a 
line energy at the edge. In order to eliminate this edge energy, lipid bilayers 
commonly tend to form closed structures such as spheroids [7]. Nevertheless, 
they can transiently open due to various stimuli such as mechanical stresses 
and thermal instabilities. The formation of these transient pores is essential 
for regulation of PH, transmembrane electrochemical potential, and concen¬ 
trations of different molecules in the cell [5]. Additionally, transient open 
membranes are formed during electro-formation [8]. More recently, stabiliz¬ 
ing pores and control over their size have been pursued by means of electric 
fields [9], sonication [ID] , and use of edge-active chemical agents m The 
rapid progress in these techniques has attracted increasing attention to the 
study of the open lipid bilayers, including molecular dynamic simulations, as 
well as continuum mechanical treatment and numerical investigations of the 
equilibrium configurations [12113]. 

Theoretical studies of the equilibrium and stability of pored membranes 
have mainly relied on constitutive assumptions for the edge, which neglect its 
flexural and torsional elasticity. For instance, Boat and Rao [14] . Capovilla 
et al. ra. and Tu and On-Yang pH EJ considered the edge energy of an 
open lipid bilayer as a given constant. Tu and Ou-Yang [18] considered 
dependence of the edge energy on its geometry, namely geodesic and normal 
curvatures. Nevertheless, their assumptions on the form of the line energy 


2 



have not been precisely justified. 

May [T9] obtained the line energy of a lipid bilayer edge through optimiza¬ 
tion of the lipid packing at the vicinity of the edge. He modeled the edge as 
a semicylindrical micelle, and took the free energy per molecule to depend 
upon the chain length of the molecules, their cross-sectional area, and the 
strength of the interactions of the molecules with each other and with the 
surrounding solution. Although successful in obtaining the line tension, that 
framework did not capture the bending and torsional energetics of the edge. 
The gap in the literature to successfully relate the macro-scale edge energy 
to its microstructure has motivated the current study. 

The interactions between the constituent molecules of a material may be 
used to obtain the free-energy density function of that material. For instance, 
Keller and Merchant m have employed such a microphysical approach to 
extract the internal energy, surface tension, and bending energy of a liquid 
surface and to relate its bending rigidity to the molecular density and interac¬ 
tion potential. In a recent application of the work of Keller and Merchant [2tTj, 
the Canham-Helfrich free-energy density for a lipid vesicle was derived based 
on microphysical considerations m Using the same approach, a model for 
the elastic free-energy of wormlike micelles was derived [[22]. In so doing, 
the surfactant molecules comprising the wormlike micelle were assumed to 
have constant length, and thus, were modeled by one-dimensional rigid rods. 
The resulted expression for the free energy was found to be quadratic in the 
curvature and torsion of the centerline of the micelle m 

The current study adopts the microphysical approach of Keller and Mer¬ 
chant (20| to investigate the elastic behaviour of the edge of a lipid bilayer. 
Following May jTU] and motivated by previous studies [23], 121, 2T>> [2H], the 
edge is modeled as a semicylindrical surface. In addition, the phospholipid 
molecules comprising the edge are modeled as one-dimensional rigid rods of 
constant length, oriented perpendicular to the centerline of the edge. The 
applied framework enables us to extract the form of the free energy and 
the flexural and torsional moduli of the edge, based on the intermolecular 
energetic interaction between phospholipid molecules. 

To find the free-energy density of the edge at a position x, we account 
for the interactions between all phospholipid molecules on the edge within a 
cutoff distance <5 from the molecules at x. We assume that the phospholipid 
molecules are perpendicular to the centerline of the edge. Our derivation 
relies on Taylor series expansions with respect to a dimensionless parameter 
g := <5/£ <C 1, where £ is a characteristic size parameter of the edge, such 
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Figure 1: Mathematical identification of an open lipid bilayer as an open surface S with 
boundary C = dS on which a Darboux frame has been shown. Also the schematic ar¬ 
rangements of phospholipid molecules in an interior point on S and at the vicinity of the 
edge C are depicted at a point. 

as its length. For l taken as the length of the edge (or equivalently, the 
perimeter of a pore), it can be related to the thickness or the length of the 
constituent molecules, if the density of the molecules along the edge and 
their aspect ratios are provided. The net free-energy of the edge results from 
integrating the free-energy density </> over the centerline of the edge. 

The paper is structured as follows. In Section [2j required mathematical 
definitions are presented. Modeling assumptions for the edge of an open lipid 
bilayer are synopsized in Section [3} Section [4] is concerned with the derivation 
of the free-energy density of such an edge. In Section |5j the consequences of 
choosing a spheroidal-particle potential (Berne and Pechukas [27] and Gay 
and Berne |28| ) are considered to obtain the material parameters present 
in the derived model. As an illustrative example, a simplified model for a 
pore on a lipid bilayer is given in Section [6j and the parameters obtained 
in Section [5] are used to find the free-energy of the pore as a function of its 
size. Finally, the key findings of the study are summarized and discussed in 
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Section [TJ Details of the various derivations are provided in the Appendix. 


2. Differential geometry of the bounding curve 
of a surface 


Consider a smooth, orientable, open surface S representing the open lipid 
bilayer, with boundary C = OS, as depicted schematically in Figure [lj Let 

C = {x : x = x(s), 0 < s < L}, (1) 

denote the arclength parametrization of the closed boundary curve C. On 
denoting the differentiation with respect to the arclength s by a superposed 
dot, it follows that cc = 1, and thus, 

x ■ x = 0, and \x x x\ — |cc|. (2) 

The unit tangent of C is introduced, in terms of the arclength parametrization 
x, by 

t := x. (3) 

Since the unit tangent t has a constant length, its arclength derivative 
t = dt/ds is perpendicular to it, and thus, perpendicular to the curve C. 
The orientation of t is called the unit normal of C, and is denoted by N. 
The curvature vector n at any point of C is then defined by the arclength 
derivative of the unit tangent t as 

k := t = fiN, (4) 

where k denotes the magnitude of the curvature of C at that point, which is 
given in terms of the arclength parametrization x, by 


K = \x X x\ = \x 


(5) 


For an arbitrary point on curve C at which n ^ 0, the unit binomial vector is 
defined by B = t x N. The unit tangent t, unit normal N, and unit binomial 
B at each point of C, form the Frenet frame {t, N, B} at that point. 

The torsion r of C is defined by B = — tN, and is expressed in terms of 
the arclength parametrization x as 


X ■ [x X x) 



( 6 ) 
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Figure 2: (a) The schematic of a section of the edge of an open bilayer, (b) cross-sections 
of the edge at positions x(s) and x(t) with Darboux frame {t,n,p} at those positions. 


The torsion r of C, describes the tendency of the curve C to move out of its 
osculating plane at a given point, or, equivalently, it measures the turnaround 
of the unit binormal B of C at a given point. In general, a space curve is 
determined up to a rigid translation, by its two locally invariant quantities: 
the curvature k and torsion r, both in terms of the arclcngth parameter s. 

On the boundary curve C of the surface S, the unit normal to the surface 
is denoted by n. Also, since x is a unimodular vector, its arclcngth derivative 
x is perpendicular to x, and thus, can be considered as a linear combination 

x = K n n + K g n x i, (7) 

of the unit normal n, and the product n x x. Notice that the unit normal n 
to the surface S is different from the unit normal N of the curve C. Further, 
let p denote the unit vector in the tangent plane of S perpendicular to the 
unit tangent t while pointing outward. We call p the unit tangent-normal. 
The set of unit normal n to the surface S at C, unit tangent-normal p , and 
the unit tangent t , which is t = n x p, form the oriented basis {t,n,p} on 
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C, known as the Darboux frame (Figure [I]). Considering that the normal N 
and binormal B of the Frenet frame of C are also perpendicular to t, they 
both lie in the plane spanned by the normal n and tangent-normal p of the 
Darboux frame. Therefore, they are related to n and p by 

N = (cos 0)n — (sin0)p, 

B = (sin 0)n + (cos 0)p, 



where 0 denotes the angle between the unit normal N of the curve C and the 
unit normal n to the surface S. Following the proof provided in Appendix 
[A] derivatives of t , n, and p with respect to the arclength s along C are 
expressed as 


' t " 


0 

K COS Ip 

—ft sin 0 


' t ' 


h 

= 

— hi COS 0 

0 

T + 0 


n 

(9) 

. P . 


k sin 0 

—r — 0 

0 


. P . 



The quantity 

r g = t + 0 , ( 10 ) 

is called geodesic torsion of the curve C on S. This quantity describes the 
rate of the rotation of the tangent plane of the surface S about the unit 
tangent to the curve C with respect to the arc length s [29]. Also, r g can be 
expressed alternatively as 


T g = x • (n x h). (11) 

Further, the curvature vector k of the curve C on the surface S is the sum of 
the normal curvature vector K n , and the tangential (or geodesic) curvature 
vector K g , i.e. 

K = K n + Kg. (12) 

The normal curvature vector K n is the projection of the curvature vector k 
along the normal n of the surface S. The geodesic curvature K g is perpen¬ 
dicular to the unit normal n to the surface, and, thus, lies in the tangent 
plane of the surface S. Hence, 


K n — (k ■ n) n, and k 9 = n x (k x to). (13) 

According to (JdJ), the curvature vector k of C is k = kN. Since 0 denotes 
the angle between N and n, the magnitude K n of the normal curvature K n 

is 

K n = k ■ n = kN ■ n = k cos0. (14) 
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The magnitude k 9 of the tangential (or geodesic) curvature vector k s , is a 
bending invariant, and is given by 


K g = — p ■ t = —np ■ N — {n x t) ■ t. 

According to t — kN, and 

n = (cos , 0)N + (sin^jB, 


the right-hand side of (15) results 


Kg = Kg (n X t), 


where k„ = Ksin^. 


(15) 


(16) 


(17) 


As mentioned earlier, the geodesic curvature vector K g at any point of the 
curve C on the surface S is the vectorial projection of the curvature vector 
k of the curve C into the tangent plane of the surface S at that point. This 
quantity is an intrinsic property of the surface, which reflects the deviation of 
the curve C from a geodesic on the surface S [30]. In general, for a geodesic, 
the geodesic curvature K g at any point is zero. Further, for a geodesic, the 
unit normal N of the curve C coincides with the unit normal n of the surface 
5, or, equivalently, the osculating plane of C at each point is perpendicular 
to the tangent plane of the surface S at that point |29j. This means that 
the Darboux frame and the Frenet frame for a geodesic are the same at any 
point. 

According to (10), and the right-hand sides of (14) and (0, the arclcngth 
derivatives of {t, n,p} in (|9]) take the form 


t K n Tl KgP , 

h = ~K n t + TgP , 

P = Kgt — TgTl. 


(18) 


3. Modeling assumptions 

The phospholipid molecules comprising the edge are allocated so that their 
hydrophilic parts lie on a thin semicylindrical surface as shown in Figure [2] 
to form a core shielding the hydrophobic tails from the surrounding solution. 
The centerline of the edge is denoted by a boundary curve C. The following 
assumptions, which are based on the previously reported observations [231 
121 ESI ED E2], are considered to model the edge of an open lipid bilayer: 





(i) ■ The phospholipid molecules comprising the edge are modeled as one¬ 

dimensional rigid rods of the same length a. 

(ii) ■ The lipid molecules are assumed to be perpendicular to the centerline 

C, residing in the plane spanned by the unit normal n and the unit 
tangent-normal p, as depicted in Figure [2] This assumption is valid 
as long as the concentration of the lipid molecules on the edge C is 
sufficiently high. 

(iii) ■ The phospholipid molecules at any cross-section of the edge have uni¬ 

form angular distribution. 

(iv) • The distribution of the phospholipid molecules at any point along C 

is denoted by the molecular density function fl > 0. In contrast to 
the angular distribution, which is assumed to be uniform because of 
symmetry considerations, the molecular distribution along C may be 
nonuniform as a result of localized curvature. 

Consider a lipid molecule at the position corresponding to s on C with 
orientation 6 measured counterclockwise from the corresponding tangent- 
normal p(s), as depicted schematically in Figure^. Let the director d(s, 6) 
denote the orientation of this molecule. By the second assumption, such a 
director can be expressed as a linear combination 

d(s,0) = (cos 0)p(s) + (sin Q)n(s), (19) 

where p(s) and n(s) denote the tangent-normal and the unit normal (to S) 
at the position corresponding to s. 


4. Derivation of the free-energy density 

In this section, the free-energy density of the edge of an open lipid bilayer is 
derived taking into account the interactions between the molecules compris¬ 
ing the edge. To do so, a microphysical approach is applied, guided by the 
work of Keller and Merchant pOj. 

Consider two molecules, with directors d and d', located respectively at 
positions x and x’ interior to C. Let the interaction energy (encompassing 
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steric, electrostatic, and other relevant effects) between the molecules under 
consideration be denoted by 

tt(x,x',d,d'). ( 20 ) 

Following Keller and Merchant [20], we assume that the interaction energy 
between two molecules separated by more than a fixed cutoff distance 5 
vanishes, in which case 

D,(x,x',d,d') = 0 if \x — x’\ >5. (21) 


In the present setting, the cutoff distance 6 is required to be small relative 
to the characteristic length £ of the edge, so that a dimensionless measure g 
of cutoff distance obeys 

£,:=)« 1 . ( 22 ) 

Hereafter, we restrict attention to interaction energies H that are of the 
form (20) but are also frame indifferent [35] . It then follows that f2(a;, x ’, d, d!) 


may depend on the positions x and x! and the directors d and d' only through 
the length \x — x'\ of the vector between x and x' , the dot products (x-x')-d 
and (x — x') d ' formed by the directors and that vector, and the dot product 


d ■ d! formed by the directors. Like Keller and Merchant [20], we assume that 
dependence of the interaction energy on the length of the relative position 


vector is scaled by the ratio g defined in (22) and, thus, that 
f2(a;, x\ d, d') = 2 Cl(g~ 2 r ■ r,r ■ d,r ■ d !, d ■ d'), 


(23) 


with r = x — x’ . The factor of two in the right-hand side of (23) is for 
simplifying later calculations. Notice that H depends explicitly on 6, whereas 


Q does not. Consequently, (21), (22), and (23) yield 


H(s 2 , PuP 2 ,Pz) = 0 if s > 


(24) 


where 

pi=sgr-d 1 p 2 = sgr-d', p 3 = sgd-d r , (25) 

with r being the unit vector corresponding to the intermolecular vector r. 

As a consequence of the foregoing discussion, the net free-energy (p net of 
the edge can be expressed as 

0net = ! 2 (y j w(s,6»)d6^jn(s)ds, (26) 
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where 


uj(s,9) = / / fl(x(s), x(t), d(s, 9), d(t, drj dt (27) 


— 7r 
2 


is the free energy due to the interactions between the molecule with director 
d(s,9) at x(s) with all other molecules and where a factor of one-half com¬ 
pensates for the double counting of interactions arising from integrating over 
both s and t from 0 to L. From (26), the free-energy density 0 at position 
£c(s) on C is simply 


— 7T 
2 


u(s, 9) d6U II(s). 


(28) 


The function fl denotes the density of the lipid molecules at any point of 
the curve C. Following the proof which relies on the Taylor series expansion 


of the integrand of (28) with respect to g up to the second derivative term, 


provided in Appendix B (28) becomes 


4> = k Q + kiK 2 + k 2 K 2 n + k 3 n g + k±n n + k 5 K n K g + k e r 2 


(29) 


which includes a quadratic expression in terms of n g , n n , and r„. Also, k 


and the coefficients kj in (29) are provided in Appendix B Notice that k Q is 


the standard line energy of the edge—the part which is independent of edge 
geometry—while the coefficients /c*, i = 1,2, 5,6 represent the flexural and 
torsional rigidities of the edge. Up to the second derivative term of the Taylor 


expansion considered here, the derived model (29) contains the linear terms 
of the normal curvature K n and the geodesic curvature n g of the boundary 
curve C, while it does not incorporate the linear term of the geodesic torsion 


r g of C. Considering a simplification of (29) in the form 


0 = ko + k^Kg - Kgo) 2 + k 2 {K r 


Hr. 


where 


Hg 0 


2k\ ’ 


^rj.o 


3 ) T k$K n Kg T kftTg , 


k,4 

2k 2 ' 


(30) 


(31) 


it can be inferred that the remaining coefficients k 3 and Aq are related to the 
spontaneous geodesic and normal curvatures n go and n no of C. This transpires 
that our model captures the spontaneous normal and geodesic curvatures H no 
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and K go of the edge, while an spontaneous geodesic torsion is absent from 
the free energy. Further, (29) and (30) include the coupling of the normal 
curvature K n 


and the geodesic curvature K g via the term k. 


n K a . However, the 


couplings of the geodesic torsion r„ with the normal and geodesic curvatures 


K r . 


and K g are absent. Our calculations show that the latter terms only 


appear by including higher order terms in the Taylor expansion, and hence, 
are of less significance. In addition, if the molecules have non-uniform angular 
distribution, i.e. the molecular distribution function n is allowed to depend 
upon 9 or r], the model would also include a linear term in geodesic torsion 
T g that would lead to the presence of a spontaneous geodesic torsion. 

The first term k Q on the right-hand side of (30) is insensitive to the shape 
of the boundary. Since the molecular distribution on the boundary has im¬ 
plicit dependence upon the ambient temperature and concentration, these 
effects may be encompassed in k Q and in the moduli k\-k§. The line tension 
k 0 has been obtained through experiments and molecular dynamic simula¬ 
tions for various types of lipid bilayers (see [23J for a comparison of different 
measurements). However, there is not enough literature on measurement of 
the remaining coefficients. By fitting our model to existing measurements of 
the line tension, the controlling parameters of the interaction potential can 
be evaluated, and further used to obtain the remaining coefficients in (29). 

The derived model (29) can be simplified into the previously presented 
models for the free-energy of the edge. In particular, the general form (29) 
provides a development to the theoretical investigations of open lipid bilayers 
presented by Tu and Ou-Yang ra, and Guven et al. [3T] . 


4.1. Total free-energy of the edge 


The net free-energy 0 net associated with the elasticity of the edge of the 
open lipid bilayer is simply obtained by integrating the free-energy density 


0 in (30) over the centerline C of the edge by 


J net 


0 ds. 


(32) 


5. Applying a concrete interaction potential 

The interaction potential in the model developed in the previous section was 
assumed to be a general function of four frame-indifferent arguments in terms 
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of the intermolecular vector and the orientation of phospholipid molecules. 
There are numerous concrete models for such interaction potentials between 
axisymmetric particles, which are vastly employed for numerical simulations 
of liquid crystals and other similar systems. Our derivation gives rise to 
integral representations for elastic moduli k 0 -k e of the edge. Substituting 


the general form G in (B.l) with an available interaction potential yields the 


material parameters k 0 -ke appearing in (29). 


One of the standard examples among such pair interaction potentials 
is the spheroidal-particle model proposed by Berne and Pechukas [27] , and 
Gay and Berne [28], in which the molecules are approximated by ellipsoids 
of revolution, or spheroids (see Figure [3]). According to such model, the 
interaction potential between two molecules with the intermolecular vector 
r and the directors d and e possesses the multiplicative decomposition 


S)(r,d,e) =£(r,d,e)((r,d,e), 


(33) 


where £(r,d, e) and ((r,d,e) denote the strength and distance parameters 
respectively. The strength parameter £(r, d , e) depends upon the orientation 
of the molecules and that of the intermolecular vector r through (Gay and 
Berne [28] ) 


€(r,d,e) = 


4£o 


(1 — \ 2 (d ■ e) 2 )G 2 


1 - — 
2 


[r-d + r-e ) 2 (r d-r • e) 2 
1 - x'd ■ e 


1 + x'd ■ e 


(34) 


with u, //, and £ 0 the fitting parameters to be chosen. More specifically, v 
depends upon the arrangement type of the molecules (e.g. side-to-side or end- 
to-end), whereas £ 0 is a constant that specifies the kind of molecules under 


consideration. The parameter x i n (34) is the shape anisotropy parameter, 


given in terms of the ratio of the length a e to the breadth a s of the molecules 
by 

* = (35) 
(<T e /(7 s J 1 


Also, the parameter x' bi (34) is given by 

(■ Se/Ss) 1 ^ 


X = 


(e e /£ s y/» + V 


(36) 
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Figure 3: (a) The schematic of a phospholipid molecule modelled as an ellipsoidal particle. 


where £ e and £ s denote the strength parameters for end-to-end and side- 
to-side arrangement of the molecules respectively. The distance parameter 
C( r,d,e ) in (33) is given by (Berne and Pechukas [27]) 

|_|2 


CO ,d,e) = exp 


O0,d,e) 


(37) 


where s(r,d,e) is called the range parameter and is given as a function of 
the orientation of the molecules and that of the intermolecular vector r by 


s(r,d,e) = cr 0 


1-* 

2 


[r-d + r -e) 2 (r d-r e 


> 2 X =1 
2 


(38) 


1 + xd-e 1-xd-e 

In (38), a Q is related to the breadth of the molecules, cr s , via cr 0 = V2a s . 


Following Whitehead et al. [34], the parameters /j and v are chosen as 

v = —1, and n — 2. (39) 


By applying the interaction potential (33) in (B.7), and assuming a con¬ 
stant molecular density II along the boundary C, the coefficients k t in (29) 
are obtained as 


/uo 

nOoCTo 

k\ 

n 2 ^ 

+ 

k 3 = k 4 


= v^erf (x)I 0 , 
k 2 

= 


xe 


= erf 0) ( x 2j + (x 2 + 2 ) J ) 

(x 2 (J + 7)(2x 2 — 3) — 21 (2x 2 + 3) j, 


192 

= k 5 = k 6 = 0 


(40) 
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Figure 4: (a). Schematic of the constant line energy /c 0 /n 2 £ 0 Oo in terms of the dimension¬ 
less cut-off distance 5/a 0 . (b) Schematic of the flexural rigidity /ci/n 2 £ 0 CTo in terms of the 
dimensionless cut-off distance 5/a 0 . As is evident from the plots, the change in k 0 /H 2 ^ 0 a 0 
and fc 1( /n 2 £ 0 <To is negligible after S = 3a 0 . Consequently, the effective cut-off distance 
after which the potential decays rapidly, can be reasonably approximated by S = 3cr 0 - 
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where x is a dimensionless parameter defined as the ratio of the cut-off dis¬ 
tance 5 to c 0 as 


(41) 


and J 0 , / and J are integral representations shown in Appendix C Hence, 
the free-energy density (j) in (29) specializes to 


</> = k Q + + K 2 n ) = k Q + ki k 2 . 


(42) 


A single phospholipid molecule can be envisioned as a molecule in which 
a water-soluble spherical head is attached to a pair of water-insoluble tails. 
Here we rely on the dimensions of a specific kind of phospholipid molecule 
(DPPC/Water system) reported by Mashaghi et al. [3Sj According to their 
estimation, the length of the aforementioned phospholipid molecule from the 
center of the head-group to the tail is ~ 22.5-30 A°, and the diameter of the 
head-group is ~ 7-10 A°. The total volume of that molecule is thus the sum 
of the volume of the spherical head and that of a cylindrical tail-group. Based 
on the equality of the volumes of the phospholipid molecule and that of the 
spheroidal replacement, the aspect ratio cr e /a s of the spheroid in Figure [3] is 
obtained between 3 and 4. Hence, the schematic of the constant line energy k Q 
and that of the flexural rigidity k\ are depicted in terms of the dimensionless 
cutoff distance 8/a 0 in Figure [I] for molecular aspect ratio cr e /a s between 
1 and 5. According to Figure [4j the change in the constant line energy k Q 
and flexural rigidity k\ is negligible after some value of the dimensionless 
cutoff distance 8/a 0 . Therefore, a rather conservative choice for the cut-off 
distance, which guarantees inclusion of all significant molecular interactions, 
is 

8* = 3 a 0 . (43) 

This result has been used to obtain the constant part of the line energy k Q 
and flexural rigidity k\ of the edge in terms of the molecular aspect ratio, as 
depicted in Figure [5} 


6. Illustrative example: dependence of free- 
energy on the pore size 

In order to estimate the change of energy of a pore with its size, consider the 
simple case of a spherical lipid bilayer with radius R , with a pore of radius r 
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Figure 5: (a). Schematic of the constant line energy fc 0 /n 2 ^ 0 cr 0 in terms of the aspect ratio 
<j e /cr s for a cut-off distance 5/a 0 =3. (b) Schematic of the flexural rigidity k\ /n 2 £ 0 cr 0 in 
terms of the aspect ratio cr e /cr s for the cut-off distance 5 = 3 a Q . 
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Figure 6: Schematic of a pore on a spheroidal lipid bilayer 


at a distance h from its centre, as depicted in Figure [6j For such a pore, the 
total curvature of the boundary curve C is 1/r, and the geodesic torsion r g 
vanishes. Also, the normal and geodesic curvatures find the frorns 


Kn = and K g = y/ k 2 - K 2 n = 


VR 2 ~ r 2 


rR 


(44) 


As a result, the free-energy density (30) specializes to 


= k Q + ki 


Vr 2 


rR 


K 


9 ° 


+ k 2 ( — - H r 


+ k§ 


Vr 2 


rR 2 


(45) 


which with (32), yields a representation for the net free-energy </> net of the 
pore as 


y net 


2nr =k ° + kl i 


VlR 


rR 


K, 


9 ° 


2 , f l 

+ h \R 


Hr, 


+ k* 


Vr 2 


rR 2 


(46) 


It was demonstrated in the previous section that when the spheroidal interac¬ 
tion potential m is employed in the microphysical model, the spontaneous 

, and the coefficient h 5 vanish, while the bending mod- 


curvatures K go and k 


uli k\ and k 2 find the same value. Using those results in (46) yields 


0net , i^l 

— = k a r + —, 
27T r 


(47) 


whereby the energy of the pore does not depend on the size of the lipid bilayer, 
nor on the placing of the pore on it. For the cut-off distance 6* = 3cr 0 and 
the molecular aspect ratios (T e /a s = 3 and a e /a s = 4, the dependence of the 
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Figure 7: Schematic of the net free-energy </>net/n 2 £ 0 cr 0 versus the scaled pore-size r/cr 0 
for two values of the aspect ratio cr e /ij s and for the scaled cutoff distance S = 3 a 0 - 


net free-energy to pore size has been demonstrated in Figure [7} The second 
term in the right-hand side of (47) leads to a minimum point for the free 
energy for r < cr 0 , which does not fall in the physically-relevant ranges of the 
pore size. For reasonable values of r/a 0 , the first term on the right of (47) is 
dominant and the dependence of the net energy on the pore size is effectively 
linear. 


7. Discussion and Summary 

An expression for the free-energy density of the edge of an open lipid bi¬ 
layer was derived taking into account the interaction between the constituent 
molecules. The resulting expression contains quadratic terms in geodesic cur¬ 
vature, normal curvature, and geodesic torsion of the boundary curve and 
a term including the multiplication of geodesic and normal curvatures. The 
derived free-energy of the edge is the evidence of an excess energy due to 
the specific arrangement of the phospholipid molecules in the vicinity of the 
boundary of an open bilayer, in accord with the results of the existing molec¬ 
ular dynamic simulations [231 24] [25], 2(i . Further, our study supplements 
the previous molecular dynamic simulations [23l [241 25] [26] and theories [19] 
in which the free energy was obtained as a constant, by providing the con- 
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tribution due to bending and torsional energies. For certain classes of lipid 
bilayers, the bending free-energy, which can be captured by our framework, 
is of more importance in contrast to the line energy [361 unj. In addition, 
our microphysical model justifies the constitutive assumptions that appear 
in continuum mechanical theories for open lipid bilayers [121 [35] [HTj . 

Our derivation gives rise to integral representations for the material pa¬ 
rameters present in the model. Specifically, the molecular origins of the 
spontaneous curvatures of the edge of a lipid bilayer have been investi¬ 
gated. A concrete soft-core interaction potential for axisymmetric rod-like 
molecules was applied on the derived model to obtain those material pa¬ 
rameters. Hence, a special form of the interaction potential suggested by 
Berne and Pechukas [27] was employed to further explore the developed mi¬ 
crophysical model. Assuming that the molecules are uniformly distributed 
along the edge, (i.e. n ^constant), the spontaneous curvatures and the tor¬ 
sional contribution to the energy vanish, resulting in (42), which includes 
only a constant part k Q , which can be interpreted as a line tension, and a 
contribution due to bending with flexural rigidity k\ = k- 2 - The dependence 
of the parameters k 0 and k± on the aspect ratio a e /a s and the dimension¬ 
less cutoff distance <5/cr Q , for various aspect ratios common for phospholipid 
molecules comprising lipid membranes, was investigated. It was concluded 
that increasing the cut-off distance after a value 8* = 3 a 0 does not affect 
those parameters. This result and the definition <j 0 = \f2a s reveals that each 
molecule on the edge interacts with less than 6 molecules in its vicinity along 
the edge. In view of this observation, the cut-off distance was set to 8* to 
study the dependence of the line tension k a and flexural rigidity on the 
molecular aspect ratio a e ja s . It is evident from Figure |4] that the line tension 
k Q is not as sensitive as the flexural rigidity k\ to the molecular aspect ratio 
(T e /cr s . This difference can be quantified by considering the relative change 
of those quantities associated with a same increment of the molecular aspect 
ratio a e /a s . Further, with an increase in the aspect ratio, the line tension 
decreases while the flexural rigidity increases. This signals that for molecules 
with greater aspect ratios, inclusion of the bending contribution to the free- 
energy is of more significance. Considering that a rod with a larger diameter 
shows more resistance to bending, the observation that the flexural rigidity 
of the edge is greater for larger molecular aspect ratios (longer molecules lead 
to larger cross-sectional diameter of the edge), agrees with what would be 
expected intuitively. 

The energy functional obtained in (30) was used to explore energy of 
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the degenerate case of a circular pore on a spheroidal lipid bilayer, resulting 
in the expression (47). The pore size needs to be greater than the mean 
distance of molecules cr 0 . For such sizes, the contribution of the flexural part 


to the energy given by (47) is negligible, and effectively, the energy increases 


linearly with the pore size. Inspired by earlier investigations on rupturing of 
soap films [3D], Litster [ID] developed a continuum model for the free-energy 
0 needed for opening-up of a pore in a lipid membrane in the form 


Y -p 1 2 

— = i r - 7 s r , 

2vr 2 


(48) 


with T the line energy of the edge and the interfacial surface tension. It can 


be inferred from (48) that transient pores with sizes less than a critical radius 
r* = T/ 7 S tend to reseal, while those having the size exceeding this critical 


radius, might grow indefinitely, leading to rupture of the membrane 
In other words, the energy to form a pore of radius r is determined by 
a balance between two competitive contributions: the energy required to 
create the edge of the pore, and the energy released by the pore surface [42]. 
Nevertheless, it is a common knowledge today (see for example the review 
by Jahnig 031) that lipid membranes possess zero surface tension, by which 
the second term on the right hand of (48) vanishes. This transpires the 


agreement of the current result based on molecular interactions, with that 
obtained previously on continuum grounds. The increasing cost of generating 
a larger pore confirms the stability of a lipid membrane with respect to the 
fluctuations that might bring about transient pores. Furthermore, the growth 
of stable pores in homogenous lipid bilayers is only possible in the presence 
of external stimuli such as an electric held. 


The class of the interaction potentials (23) selected in the present study 


and the integration (26), only account for the lipid bilayers in which the 


physiochemical properties of the constituent molecules are identical. An im¬ 
portant corollary of our model would follow from a generalization of the 


arguments of the interaction potential (23) and the integration (26), to allow 


for the interactions between phospholipid molecules of different physiochem¬ 
ical properties. Such a generalization permits modeling perforated mixed 
lipid bilayers, such as those reported by Ogl§cka et al. m and Jiang & 
Kindt [36] . Further, our model accounts for the elastic free-energy of the 
edge of open lipid bilayers in which the lipid molecules are tilted only at 
the edge, forming a semicylindrical rim along it. Another generalization of 
the present work would include tilt fields of smaller gradient, such as those 
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considered by Hamm and Kozlov [32], Rangamani and Steigmann 1415], and 
Rangamani et al. m- In such an approach, the gradual tilt at the vicin¬ 
ity of the edge changes the thickness of the lipid bilayer and, thus, leads to 
the deviation of the conformation of the edge from a semicylindrical shape. 
These potential generalizations remain to be investigated in future. 
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Appendix 


A. 


Detailed derivation of 



Let 0 denote the angle between the unit normal N to the curve C and the 
unit normal n to the surface S. The Darboux frame {t. n. p) at any point on 
C, is obtained by rotating the Frenet frame (t, N, B} about the unit tangent 
t by that angle. Hence, 


' t ' 


n 

= 

. P . 


t n 

i T 

P\ 


0 0 

cos 0 sin 0 
— sin 0 cos 0 



t 


N 


B 


(A.l) 


denote the Frenet frame at a given point on C, where T denotes the transpose. 
Also let 

'10 0 

Q := 0 cos ip sin 0 (A.2) 

0 — sin 0 cos 0 

denote the transformation between the two frames. Thus, 


X = QY, or Y = Q t X. 


(A.3) 


The derivative of the Frenet frame (t, N, B} of C with respect to the arclcngth 
s, follows from the Frenet-Serret formulas [3E] 


t 


0 K 0 


t 

N 

= 

1 

o 


N 

B 


0 -r 0 


B 


(A.4) 
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Thus, 


Y = AY = AQ t X, 


where 


A = 


0 k 0 
—k 0 r 
0 -r 0 


Taking arclength derivative from both sides of (A.3)i yields 


X = QY + QY, 


(A.5) 
(A.6) 


(A.7) 


where 


Q 


"0 0 0 

0 —0 sin 0 0 cos 0 

0 —0cos0 — 0sin0 


(A.8) 


In view of (A.30 and the right-hand side of (A.5), (A.7) takes the form 


X = (QQ t + QAQ T ) X, 


(A.9) 


or, equivalently, 


' t ' 


0 

K COS Ip 

—k sin 0 


' t ' 


h 

= 

— K COS 0 

0 

r + 0 


n 

(A. 10) 

. P . 


k sin 0 

—r — 0 

0 


. P . 



In view of K n = kcos0, K g = n sin 0, and T g = r + 0, (A. 10) simplifies 
to (18). 


Appendix B. Derivation of the free-energy den¬ 


sity (29) 


In this Appendix, the expansion of (28) to (29) is presented. As mentioned 


earlier, only the molecules separated by a distance less than 5 may interact. 


Hence, the domain of the integral with respect to t in (27) is replaced by 
[—5,5]. Upon replacing uj(s,9) in (28), and substituting H by its equivalent 
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from (23), and applying the change of variable t — t Q = sg , (28) takes the 
form 


= Q 



! —P.J ~ 7r ,/ ~ 7r 
2 2 


g r r,r ■ d(t 0 , 0 ), 

r ■ d(t D + sp, rj),d(t 0 , 6) ■ d(t D + sp, p) j 

n(t 0 )II(t 0 + sp) dridOds. (B.l) 


It is necessary to expand the right-hand side of (B.l) in powers of p neglecting 
terms of o(p 2 ). Introducing the abbreviations 


n:=n( 0), t := t( 0), p := p(t 0 ),' 

n : =n(o), ri:=n(o), ii:=n(o), 


(B.2) 


and applying the identities given in (18), the following expansions up to p 2 
are obtained: 

x{sg ) = £c(0) + gst + ^ s 2 g 2 K n n - ^s 2 g 2 K g p + o(p 4 ), 

n(sg) = (^-(k n - K g Tg) — sgKy^j t + (l ~^-( K l + T g) S jn 

+ (^-{KnKg + fg) + sgr g ^jp + o(g 2 ), 

fs 2 g 2 \ ( s 2 g- \ 

P(sg ) = f — (hi n K g - r g ) - sgr g j n + ( gsn g + (k 9 + T g K n ) j t 

+ (l-^(r 2 + K 2 g )^jp + °(g 2 ). (B.3) 


Therefore, the arguments of the interaction potential in the right-hand 


side of (B.l) become 

g~ 2 \x(t 0 ) - a?(t e )| 2 = s 2 + A 1 g 2 s 4 + o(p 2 ), 

(x(t 0 ) - x(t e )) ■ d(t 0 ,6) = A 2 g 2 s 2 + o(p 2 ), 
(x(t 0 ) - x(t Q )) ■ d(t e , rj) = A 3 g 2 s 2 + o(p 2 ), 
d(t a , 6) ■ d(t e , rj) = A 4 + A 5 gs + A 6 g 2 s 2 + o(p 2 ), 


(B.4) 
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where t e = t a + sg, and 

k 2 a K g cos ® ~ K n sin 9 K n sin rj — K g cos 77 

Al —-, Ao — -, Art == -. 

12 ’ 2 ’ 2 

A 4 = 008(6* — rj) A 5 = — T g sin(6* — 77), 


K n Kg 

Ae, = 


T T 

^ sin(# + 77 ) —^ sin(d — 77 ) —cos (0 — 77 ) 


k 2 K 2 

fl 'g n n ■ n ■ 

-- cos a cos 77-sm t) sin 77. 

2 2 


Expanding £2 up to g 2 and using 




n(s^) — n + s^n 4———n + o(g 2 ), 

results in the following energy-density for the edge 


(B.5) 


(B. 6 ) 



£>£2(s,#,77)n 2 d77d# ds + -g 3 Q(s,6,r))s 2 UIl > dr) &6 ds 


£ 2 i(s, 9, g)s 2 + 6 £ 2 4 (s, 9 , rj) cos 9 cos 77 ) dr) c!6* ds 


£2i(s, 9 , 77 )s 2 + 6 fl 4 (s, 0, 77 ) sin 9 sin 77 ) d ?7 d 9 ds 


£2 2 (s, 9, rj) cos 9 — £2 3 (s, 9 , 77 ) cos 77 ) d ?7 c! 6 * ds 


fl 3 (s, 6 *, 77 ) sin 77 — 0 2 (s, 6 *, 77 ) sin 0 ) d /7 c! 6 * ds 


£ 2 44 (s, 9 , rj) sin 2 (d — rj) — fi 4 (s, 9, rj) cos (9 — 77)^ dr) d 9 


£ 2 4 (s, 9 , rj) sin(6* + rj) dr) d9 ds >, 


(B.7) 


£ 2 (s, 9, rj) := 0 (s 2 , 0 , 0 , cos (9 — 77)), 
£2j(s, 9, rj) \= Cl ti (s 2 , 0, 0, cos(d - 77)), 
&u(s,9, 77 ) := Cl,a (s 2 , 0, 0, cos(6* - 77)),. 


7 €{1,2,3,4}. 


(B. 8 ) 


ds 
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Note that fl, J2, and 0 M vanish when s > l. Thus, the final expression for 
the free-energy density of the edge becomes 


(j) = ko + kiK 2 g + k 2 K 2 n + k 3 K,g + k±K n + k b K n K g + k§T 2 , (B.9) 


or, equivalently, 

(f) k 0 T k\(hvg + ^g(^n ^no) d~ k^Hv n hvg + k^Tg , 

where the parameters k 0 -k§ are given by the integral representations in 


(B.10) 

Eg. 


Appendix C. Integral representations / 0 , / and 

J in (| 40 |): 


Io 
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X 2 cos 2 (6 — rj) drj dd, 




1 

y/l — X 2 COS 2 (6 — 7 j) 

cos(2 6 — 2 rf) 
y/l — X 2 cos 2 (9 — 7]) 


d'l] d 9, 
d ■/? d 9. 
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